Analytical expressions for the polarizability of the honeycomb lattice 
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We present analytical expressions for the polarizability P fl (g :E ,aj) of graphene modeled by the 
hexagonal tight-binding model for small wave number q x , but arbitrary chemical potential /j,. Gen- 
erally, we find P l _ l (q x ,uj) — P^(u>/uj q ) + q x P^"(uj) with u) q — v F q x the Dirac energy, where the first 
term is due to intra-band and the second due to inter-band transitions. Explicitly, we derive the 
analytical expression for the imaginary part of the polarizability including intra-band contributions 
and recover the result obtained from the Dirac cone approximation for /x — > 0. For /i < V3t, there 
is a square-root singularity at uj q — v F q x independent of /j. For doping levels close to the van Hove 
singularity, fi = t ± 8fi, ImP M (q x , ui) is constant for 8[i/t < uj/uj q <C 1. 



PACS numbers: 81.05.ue, 73.22.Pr, 74.70.Wz 



I. INTRODUCTION 

Graphene is a two-dimensional carbon allotrope which 
has attracted immense research activity due to its novel 
mechanical and electronic properties. 1-4 It is also inter- 
esting in view of potential applications in nanomechanical 
and nanoelectronic devices and one of the reasons for this 
lies in the possibility to change the carrier density and 
type by applying a gate voltage between graphene and 
the isolating substrate. This gave rise to the celebrated 
ambipolar field effect with an almost constant mobility 
of around 10 4 cm 2 /Vs for graphene on Si02-substrate. 5 

But it is also possible to modify the Fermi surface via 
chemical doping, i.e., by the deposition of e.g. potassium 
atoms which donate their lone valence electrons to the 
surface layer. 6 Whereas with a typical back-gate voltage, 
chemical potentials of the order of a few 100 meV can 
be reached, doping levels up to the van Hove singularity 
were recently achieved via chemically n-type doping us- 
ing various combinations of K and Ca on both sides of 
graphene. 7 

For small energies and chemical potentials w,/i^leV, 
graphene can be well described by the Dirac cone ap- 
proximation, optionally including weak trigonal warping 
corrections. 8,9 For larger energies or doping levels around 
the van Hove singularity, the linearization around the K- 
pionts is not applicable anymore and one usually has to 
resort to numerical methods. 10 Only around special sym- 
metry points like the M-point, analytical solutions are 
still feasible. 11 - 12 

The tight-binding model on a hexagonal lattice with 
only one orbital is a good approximation to the electronic 
structure of graphene for doping levels up to the M- 
point. 13 Including interaction effects via exchange self- 
energy corrections preserve the trigonal warping of the 
Fermi surface topology 15 and only close to the van Hove 
singularity the local Coulomb interaction U can lead to 
instabilities at low energies. 14 Superconductivity medi- 
ated via electron-electron interaction was hence predicted 
for graphene provided that the Fermi surface is close to 
the M-point. 7 

Instabilities due to Coulomb interaction can be an- 



alyzed in terms of the effective interaction vertex 
which is dressed by the static polarizability as a first 
approximation. 16 Usually this is done numerically and 
only for gated graphene with small chemical potential, 
analytical expressions were obtained within the Dirac 
cone approximation. 17-21 Here, we will report on an an- 
alytical solution for the imaginary part of the polariz- 
ability for arbitrary chemical potential within the tight- 
binding model. The solution is valid for low energies 
huj<2fi and small incoming wave vector |q| = q x where 
the x-axis denotes the high-symmetry direction which 
connects the T- and the M-point. Analytical expressions 
are especially important close to the van Hove singular- 
ity since numerical methods tend to fail in this regime. 
They also provide insight in the scaling properties and 
can explain the peak-splitting in arbitrary direction of 
the square-root singularity seen in the numerical solution 
ofRef. 12 . 



II. POLARIZABILITY OF THE HONEYCOMB 
LATTICE 



The density-density correlation or Lindhard function 
of the honeycomb lattice is given by 

P^u) = -%( d 2 k *Xv(k,q) 

W Jl.BZ xx , =± 

n F (E x (k)) -n F {E x '(k + q)) 



F\,\' (k, q) = f\.\>(k,q) 



E x (k) - E x '(k + q) + hcu + iO ' 

(1) 



with the eigenenergies £ ;± (k) = ±i|</>k| (f ~ 2.7eV is the 
hopping amplitude), np(E) the Fermi function, g s = 2 
the spin-degeneracy and 0k = ^2 S e lk S the complex 
structure factor where 5 denote the three nearest neigh- 
bor vectors. Due to the two gapless bands, the above 
expression contains the band-overlap function 



/±(k,q) = i[ I ! R- 



y k+q 



k+ql 



(2) 
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FIG. 1: a) The hexagonal and rhombical Brillouin zone b) The 
symmetrized rhombical Brillouin zone and its segmentation. 
The inner square refers to j = — , the outer triangles refer to 
j = +; additionally the values of s and s' are given. 



which marks the crucial difference to standard text-book 
results containing only one band. 

It is possible to transform the integration over the 
rhombical Brillouin zone into an integration over a 
quadratic area. First, we symmetrize the Brillouin zone 
by introducing dimensionless parameters 3ak x /2 — »■ k x 
and \/3ak y /2 — > k y , with a = 1.4A the carbon-carbon 
distance. Second, the integration over the outer triangles 
of Fig. 1 is written as the integration over the quadratic 
area by symmetrizing the sign of the external wave num- 
ber q — » (sq x ,s'q y ) with s,s' = ±. The integration over 
the inner squares of Fig. 1 can be written as the inte- 
gration over the quadratic area by slightly modifying the 
integrand i*\,A'(k, q). For this we introduce the gener- 
alized energy E X ' J (k) = A|^| and overlap function f± 
with 



3 + 2 cos(2/c y ) + jA cos(fcy) cos(fc x ) (3) 



and 



/i(k,q) = - 1± 

2 \ |0k||0k+q| 



(4) 



where 



f 3 (k, q) = cos(q x ) + j2 cos(k y ) cos(k x + q x ) (5) 
+ 2 [2 cos(fc y ) cos(<7;r/2) + j cos(k x + q x /2)} cos(k y + q y ) 

and the dimensionless variables q x a — > q x , \f3q y a/2 — > 
Qy 

Changing now the integrand i*\,A' (k, q) to F 3 X x , (k, q) 
by substituting E x (k) -> E x ^(k) and /±(k,q) -> 
/±(k,q), we have 

x E E E Fl,y^sQ^s'q y ) . (6) 

A,A'=± j=± s,a'=± 

From this definition of the polarizability, we can make a 
suitable substitution to obtain an analytical solution. 



III. ANALYTICAL SOLUTION 

For finite chemical potential \i > 0, zero temperature 
T = and small wave vector |q|a <C 1, the imaginary 
part of the polarizability for low energies fiw<2/^ is only 
determined by intra-band transitions. This limit shall 
be denoted by ImP^ — > ImP< With q g)S / = (sq x ,s'q v ), 
/i = n/t, and ui = Hui/t, we obtain 



ImP<(q,w) 



,V3 



Airt 



-it/2 



dk x 



IX 1 1 



dk 



> E E 

s,s' = ± j=± 



/|(k,q 5iS /)q SiS /V k |^|J (/I - |</> k |) 5 (o>- q s , s 'V k | 



(7) 



For q-vectors in rc-dircction, i.e., q y = 0, the integral 
over the two delta-functions can be performed analyti- 
cally with the substitution x = sink x , y = smk y . 

The analytical expression depends on the zeros of the 
first delta function 



- A 2 , ft 

4 2 



1 



(8) 



where oj q = vpq x is the Dirac energy with vp = ^at/h 
the Fermi velocity. The final result can then be written 



in the compact form 
9 s >/3 



ImPMg^w) 



8tt (ftv F ) 2 



(9) 



x(9(w+-w) 



9{u 



with = #(ft — 
uj*8(jl — \/3) where w 



l)w* 



y-\ P-- v- 



and w+ 



1-2/3 



■v/To" 



ft) 



ft 2 



9/ft 2 . Note that 



ImF<((j x ,w) only depends on the ratio uj/co q . 

Real and imaginary part of a response function arc 
related via the Kramers-Kronig relation. For this, the 
w-dependence for the whole spectrum is needed, but we 
can define the contribution to the real part of the polar- 
izability that originates from intraband transitions: 



2 ["« 

ReP<(q x ,u) = - dijj 
t Jo 



,o/ImP<( ga ,h/) 

(UJ 1 ) 2 -UJ 2 



(10) 



Also ReP<(q x ,u) depends only on uj/u q . For u = 0, 
there is no dependence on q x and we have 



ReP<(q x ,u = 0)=p{fi) 



(11) 



with = X)k^(M ~ t\^k\) the density of states of 
the tight-binding model which can be expressed in terms 
of the complete elliptic integral of the first kind. 22 

Let us now discuss the f-sum rule for small q x defined 



as- 
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2 h 2 A r 



7t (aq x ) 2 t 



6t 



dojujImPft (q X7 uj) =I t 



v ■ 



(12) 



with T n 



2N 



6k I) where A = A C N is the 



area of the sample (A c — 3\/3a 2 /2). 

Let us divide the total weight into partial weights 
due to intra-band and inter-band contributions 1^ = 
!<+!>. With Eq. (9) and ft < 1, X< is then given 
by 



(13) 



where = (3 — ft, ) /2ft and £_ = (1 + ft ) /2ft which can 
be expressed in terms of elliptic functions. For ft > 1, 
the integration bounds have to be modified. 

Since 1^ and Z< are independent of q x and the lead- 
ing order of the band-overlap of the inter-band transition 
is proportional to q x , we have for the inter-band con- 
tributions to the polarizibility ImP//(q x ,uj) = q x 9(hui — 

2jLt)P(w), where P only depends on the energy. Obtain- 
ing the real part via the Kramers-Kronig relation, we can 
thus generally write the polarizability as 



P^q x , U ) = P<(uj/u q ) + <£P>(w) , 



(14) 



where the first term is due to intra-band and the second 
term due to inter-band processes. 

In table I, the relative weight of the intra-band tran- 
sition is shown as they contribute to the total weight. 



ft 


0.1 


0.5 


0.9 


1.0 


1.1 


1.5 


2.0 


2.5 




1.574 


1.558 


1.462 


1.395 


1.319 


1.076 


0.755 


0.396 


2> 


0.053 


0.266 


0.528 


0.654 


0.773 


0.927 


0.981 


0.997 



TABLE I: Total and relative weight due to intra-band tran- 
sitions of the f-sum rule for various chemical potentials fj,. 



IV. DISCUSSION 
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FIG. 2: (color online): The imaginary part of the polarizabil- 
ity as function of the energy ui for various chemical potential 
ft, = \i/t divided into the three regimes: ft < 1 (left panel), 
1 < ft < v3 (middle panel), and ft > \/3 (right panel). The 
result of the Dirac cone approximation is also shown (dashed 
line). 
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FIG. 3: (color online): The imaginary part of the polarizabil- 
ity ImP (|q|, ip, to) with |q|a = 0.01 as function of the energy 
uj for various angles ip and two chemical potentials \ijt — 0.05 
(left) and \i/t = 0.5 (right) at k B T/t = 0.01. Also shown the 
result from the Dirac cone approximation at zero temperature 
(dashed line). 



A distinct signature of non-interacting 2D electrons is 
a divergent behavior of the polarizability at the thresh- 
old for the excitation of electron-hole pairs determined 
by the Fermi velocity. For small chemical potential, the 
Fermi velocity vf in graphene is independent of /x and 
the divergence takes place at cu q = Vpq x . But even in 
the regime where the Dirac cone approximation does not 
hold, i.e., curvature in form of trigonal warping has to be 
taken into account, we find the same divergent behavior 
at uj q . It is remarkable that this threshold is independent 

of the chemical potential \i up to < 
For /j, — > 0, we obtain 

o« u ojQ{uj n — lo) 



This expression is also obtained from the Dirac cone ap- 
proximation in the limit uj, q <C /i- This shows that the 
Dirac cone approximation is valid only in the limit /i — > 
which leads to the criticality of Dirac Fermions at zero 
gate voltage and zero temperature. 24 

Let us now discuss the special case where the chemical 
potential lies at the van Hove singularity. For /i = t, we 
have 

Whereas ImP At= t(g K ,aJ — > 0) cx const is a necessary 
condition to yield the logarithmic singularity of the den- 
sity of states at the van Hove singularity, we shall now 



4 



investigate the behavior for chemical potentials close to 
the van Hove singularity /i = t±5/i. Neglecting quadratic 
contributions 0(5fi 2 ) which is valid for 5jl <C (co/oJ q ), 
only the solution y + contributes to a finite imaginary 
part for low energies and we obtain 



Im^t±{ p (gx,w^ 0) 



g s V3t±(6n/2) 
A-k (hv F ) 2 



(17) 



We thus find the onset of non-Fermi liquid behavior for 
excitations in the direction of the highest symmetry for 
chemical potentials close to the van Hove singularity. For 
energies with <C (u>/uj q ) <C 1 and \i < t, also y_ 

contributes to the constant behavior and we obtain the 
limit value of Eq. (16). 

In Fig. 2, we show the results in three panels for the 
three regimes jl < 1 (left panel), 1 < p, < V3 (middle 
panel), and jl > (right panel). For jl < 1, the charac- 
teristic feature is the square- root singularity at co = oj , 
independent of the jl. For low chemical potential /i<0.3, 
the result agrees well with the Dirac cone approximation 
(dashed curve). For larger chemical potential, there is 
increasing weight for lower energies and for jl = 1 , ImP^ 
shows constant behavior for low energies. For chemical 
potentials with 1 < jl < ^/3, there is - in addition to 
the square-root singularity at u = uj q - also a square-root 
singularity atw = w*. For jl > only the square-root 
singularity at u = uj* survives. 

For arbitrary direction ip = tan -1 (q y /q x ), in particu- 
lar for q in y-direction, the solution of the zeros of the 
delta-functions of Eq. (7) involves polynomials of 4th 
and 6th order and we were not able to obtain an ana- 
lytical solution. The numerical solution, though, shows 
the expected peak splitting of the square-root singularity 
at u = ui q , as shown in Fig. 3, with the maximal peak- 
splitting for ip = n/6 + mr/3, n e N. As suggested by the 
analytical expression, the double-peak structure appears 
even for small chemical potentials, e.g., fi/t = 0.05, for 
which the Dirac cone approximation should hold. 

The curves of Fig. 3 were obtained for ksT/t = 0.01, 



thus slightly larger than room temperature. The curves 
for ip = are basically unaffected by temperature, but 
for arbitrary direction the algebraic divergences seen for 
n/t = 0.5 are smeared out at larger temperature as it is 
the case for fi/t = 0.05. The curves for n/t = 0.05 and 
<p 7^ develop the algebraic divergence for decreasing 
temperature, so generally, we can say that the algebraic 
divergences become broadened when the energy set by 
the temperature is much larger that the maximal peak- 
splitting at ip = it/6. For the curves of Fig. 2, we find 
from the numerical solution that the singularity at oj = 
uj q is practically unaffected by temperature whereas the 
singularity at lo = lo* is strongly broadened for T > 0. 



V. CONCLUSIONS 

In summary, we have presented analytical expressions 
for the imaginary part of the polarizability for wave vec- 
tors in the T ~ M direction. The results are valid for all 
chemical potentials [i for low energies fiw<2/x and small 
wave vector |q|o < 1, As a special feature, we find a 
singularity at the Dirac energy io q = vpq up to large 

jl < \/3. For the chemical potential close to the van 
Hove singularity, we find the onset of non-Fermi liquid 
behavior with ImP M=t ±5 Al oc const for Sfi/t < cu/uj q <C 1. 
For arbitrary direction ip, there is a peak splitting at uj q 
which prevails for small chemical potentials and only for 
ip = 0, the Dirac cone approximation is recovered for 
H -> 0. 
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